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ABSTRACT 

A  general  method  to  analyse  multi-layered  Frequency  Selecti\e  Surfaces  is  given 
which  includes  a  full  coupling  of  evanescent  and  decaying  modes  in  a  well  conditoned 
manner.  The  method  is  capable  of  analysing  arbitrary  lossy  dielectric  and  magnetic 
materials  used  as  supports  with  partially  conducting  (Salisbury)  surfaces  and  is  valid  for 
incident  plane  waves  with  arbitrary  polarisation  and  angles  of  incidence. 

The  solution  is  obtained  by  modal  matching  of  the  fields  in  the  apertures  of  each 
surface  in  a  manner  which  satisfies  the  boundary  conditions  using  Galerkin  testing.  One 
important  feature  is  the  manner  in  which  the  cascaded  matrices  are  put  together  to  avoid 
numerical  problems. 

Theoretical  predictions  have  been  made  on  a  wide  range  of  surfaces  and  tested 
against  experiment  with  satisfactory  results. 
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A  WELL  CONDITIONED  FORMULATION  FOR  THE  ANALYSIS  OF 
MULTI-LAYERED  FREQUENCY  SELECTIVE  SURFACES  USING 
A  CASCADED  MATRIX  APPROACH 

Andrew  Mackay  and  J  G  Gallagher 


INTRODUCTION 

Frequency  selective  surfaces  have  many  possible  uses  as  filters,  polarisers  and  for 
radomes  and  in  other  electromagnetic  structures.  In  many  instances,  single  layer  surfaces 
might  be  suitable  and  for  many  others  a  multi-layer  surface  can  be  designed  from  a 
knowledge  of  the  Transmission  and  Reflection  coefficient  of  the  single  surfaces  in  isolation. 
This  latter  approach  is  often  referred  to  as  a  simple  scattering  matrix  (SMA)  analysis  and 
is  valid  provided  there  is  no  evanescent  mode  coupling  between  layers  [1],  However, 
there  will  be  many  instances  where  evanescent  mode  coupling  must  be  taken  into  account 
when  the  layers  are  a  sufficiently  small  fraction  of  a  wavelength  apart.  Indeed,  there  is 
no  reason  to  suppose  that  multi-layer  frequency  selective  surfaces  can  not  be  explicitly 

designed  to  utilise  this  strong  coupling.  In  such  a  design  procedure  it  is  clearly  desirable 

to  have  a  method  for  calculating  the  properties  of  such  strongly  bound  layered  structures 
with  a  fair  degree  of  confidence. 

Shuley  [2]  considered  a  generalisation  of  the  SMA  to  allow  for  complete  interaction 
by  the  use  of  evanescent  mode  coupling.  Similar  methods  have  been  used  later  by  Cwik 
and  Mittra  [3],  Hall,  Mittra  and  Mitzner  [4]  and  Christodoulou  et  al  [5]. 

In  order  to  use  an  exact  Floquet  modal  analysis  it  is  necessary  to  construct  a  unit 
cell  common  to  all  layers  of  the  surface  so  that  the  entire  structure  is  periodic  with 
respect  to  this  global  unit  cell.  This  is  only  unnecessary  if  there  is  insignificant 

evanescent  mode  coupling  between  surfaces  and  in  this  case  surfaces  do  not  need  to  share 
a  common  periodicity.  In  general,  any  periodic  surface  can  be  constructed  from  a 

(global)  unit  cell  represented  as  a  parallelogram  illustrated  in  fig  1 .1  below 
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All  fields  generated  by  the  FSS  must  thus  be  periodic  with  respect  to  the  Floquet  modal 
values 
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where  p,q  are  integers  of  either  sign,  k0  =  2ir/X0  is  the  wavenumber  of  the  incident  field 
in  free  space,  and  fl0  and  c’o  are  angles  of  incidence  of  the  plane  wave  as  illustrated 
in  Fig  1.2  below: 


FSS  in  x-y  plane 


z 


It  should  be  noted  that  these  Floquet  modes  are  common  to  ah  the  layers  of  the  FSS, 
independent  of  any  dielectric  or  magnetic  material  values  since  the  value  of  is 
unchanged  within  the  structure  and  k.sintf  is  unchanged  by  Snell's  law  throughout  the 
structure  (where  k  is  the  wavenumber  in  any  layer  of  the  multilayer  structure,  and  fl  is 
the  direction  of  incidence  -  complex  for  a  lossy  material  -  within  that  layer). 

The  modal  fields  have  an  x,y  dependence  proportional  to  exp(-j[UpqX  +  Vpqy])  and 
a  z-dependence  proportional  to  exp(-j-ypqz)  where  r 
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and  when  an  implicit  time  dependence  e-M  is  assumed.  er,pr  are  possibly  complex  with 
a  negative  imaginary  part  representing  loss  and  ->pq  is  taken  so  that  the  imaginary  part 
Im(7pq)  <  0  to  ensure  correct  behaviour  as  z  -*  +w  (assuming  all  distances  are  taken  in  a 
positive  sense). 

The  electric  field  may  be  decomposed  into  orthonormal  Transverse  Electric  and 
Transverse  Magnetic  vector  Floquet  Modal  Fields  ^pqr  with  r  =  1  representing  TE  and  r 
=  2  representing  TM  modes,  so  that  f  ’ 
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We  apply  the  electric  and  magnetic  field  boundary  conditions  by  considering  convenient 
canonical  structures  with  which  we  build  the  complete  multi-layer  structure.  For  example, 
an  FSS  interface  backed  by  a  single  layer  of  material  (as  used  by  Chen  [6])  may  be  used. 
Alternatively,  only  a  single  interface  may  be  used.  The  latter  approach  is  adopted  for 
programming  convenience  which  will  be  efficient  whenever  any  dielectric  layer  is 
sandwiched  by  FSS  conducting  structures  on  both  sides.  Choosing  this  approach,  any 
arbitrary  multi-layered  FSS  may  therefore  be  constructed  from  the  two  canonical  structures 
illustrated  below. 


where  a  metal  FSS  pattern  is  bounded  on  either  side  by  materials  with  different  bulk 
electrical  characteristics  for  fr  and  pr.  (Resistive  patterns  may  also  be  considered  -  see 
Hall,  Mittra  and  Mitzner  [4]) 


where  there  is  a  resistive  (possibly  not  present)  layei  cf  zero  thickness  bounded  on  either 
side  by  two  different  electrical  materials 
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These  two  structures,  when  separated  by  any  thickness  of  uniform  material  in  any 
combination  can  be  used  to  represent  any  multi-layer  FSS  of  conducting  elements,  resistive 
screens  and  arbitrary  supporting  materials,  as  illustrated  in  the  example  figure  below,  Fig 
2.3. 


Figure  2.3 

where  each  interface  may  be  represented  as  an  example  of  a  canonical  structure  (1)  or  a 
structure  (2).  In  this  paper,  we  choose  to  propagate  the  electric  field  through  from  one 
layer  to  the  next. 

We  now  show  how  to  relate  the  (electric)  fields  on  either  side  of  the  interface  of 
a  canonical  structure. 

BOUNDARY  CONDITIONS  AT  AN  FSS  INTERFACE 

We  first  consider  the  interface  of  canonical  structure  (1).  Using  a  generalised 
scattering  matrix  approach,  the  field  to  either  side  of  the  interface  can  be  written  as  the 
sum  of  incoming  and  outgoing  modal  fields. 
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where  x  is  evaluated  at  the  interface  z  =  0  and  the  metal  pattern  can  be  resolved  into 
regions  of  conductor  (no  aperture)  and  regions  of  no  conductor  (aperture).  Similarly,  on 
the  right  of  the  interface,  the  total  tangential  electric  field,  H+l  may  be  written  as 


x  f  FSS  aperture 

x  4  FSS  aperture 


(7) 
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pmax  and  Qmax  are  sufficiently  large  and  where  Apqr,  Rpqr,  Tpqr,  Wpar  are  complex 
constants  representing  incident  reflected  transmitted  and  returning  wave  mode  coefficients 
respectively.  We  may  now  apply  the  tangential  boundary  condition  across  a 
magnetic/dielectric  interface;  Et  is  continuous  over  the  aperture.  Since  the  tangential  field 
is  zero  elsewhere  on  both  sides,  we  can  write 
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combining  (9)  and  (10)  we  now  have 


y  a  i(i)  £  (x)  -  y  r  £(i)  £  (jt> 

A  Pqr  pqr  *-pqr  -  Pqr  Pqr  Pqr  c 


y  T  { (2)<e  (24 )  -  y  w  £^2^i2  (2£)  ,  24  f  aperture  only 

L  pqrspqrwpqr  L  pqr'pqr^pqr^  3 


Equations  (8)  and  (13)  may  now  be  solved  in  the  following  manner.  Firstly,  the 
orthogonal  properties  of  ^pqr  may  be  used  in  conjunction  with  (8)  to  write 


A  +R  -T  +W  V  P,q,r 

pqr  pqr  pqr  pqr 
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(19) 


Et<s)  -  !  FrnneW*>  '  *  ‘  aperture 

mnC 

Note  that  the  inner  product  (18)  implicitly  assumes  that  Et(x)  is  zero  when  x  <  aperture, 
and  thus  (19)  requires  imn()(x)  only  to  be  defined  over  the  aperture.  In  practice,  vmnf 
may  be  chosen  as  the  electric  field  slot  modes  for  a  slotted  structure  (or  as  the  fields 
caused  by  the  expected  current  modes  in  an  element  structure)  thus  ensuring  as  few  vmnC 
as  possible  for  good  convergence  to  the  solution.  In  order  to  solve  (17)  given  a  set  of 
tmnC  may  assume  Galerkin  testing,  substitute  (19)  into  (17),  form  the  dot  product  of 
(17)  with  k*MNI_te)  and  integrate  over  the  unit  cell  (with  defined  as  zero  for 

x  <  aperture).  This  yields  the  set  of  equations 


2(w]  -  K3  {FmnC}  -  2{w} 
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In  order  to  employ  matrix  notation  the  subscrims  (or  superscripts)  MNL  (mnC)  in 
(20)-(24)  will  be  tabled  by  a  single  index  i(j)  so  that  m  {1  i}.  {YMNLmn^}  = 

{ Yj ; }  etc.  We  may  thus  write  the  vector  of  values  1  =  {lj},  etc,  and  the  matrix  of  values 
Y  =  { Yjj },  etc.  (20)— (24 )  now  become 
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C^E^'^A  ,  E^  -  diagonal  matrix,  diag 
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cV2)W  ,  W  -  jw  }  ,  F 

so  that  (20)  may  be  rewritten  in  matrix  notation  as 

2CtE(PA  -  Cf[E(1)  +  E(2)]  C  F  -  2CtE(2)W 

Note  that  there  is  a  very  important  feature  in  this  equation;  C  is  in  general  rectangular  so 
that  neither  its  inverse  or  the  inverse  of  exist.  Indeed,  even  when  C  is  square  it  is 
likely  to  be  extremely  ill  conditioned  if  sufficiently  many  Floquet  modes  are  taken  into 
account  (ie  Pmax  and  Qmax  are  sufficiently  large).  This  seems  to  be  a  general  feature 
of  most  frequency  selective  surfaces;  the  excited  modal  fields  and  currents  always  span  a 
significantly  smaller  set  than  the  set  of  permissible  Floquet  modes  might  allow. 

Substituting  (19)  into  (16)  and  using  matrix  notation  with  {Tpqr}  ■  {T j }  ■  T,  yields 


-  W 

(25.) 


T  +  W  -  C  F  (26) 

substituting  (26)  into  (25)  gives  the  matrix  equation 

2CfE(1A  -  C  +(E(1)  ♦  E(2)]l+  cV”  -  E{2)]  S  (27) 

which  together  with  (14)  in  vector  notation, 


A  +  R  -  I  +  W 


(28) 


gives  us  the  general  connection  matrix  between  electric  fields  on  the  left  and  right  hand 
side  of  the  interface, 
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where  I  is  the  identity  matrix  and 


C 


(30) 


Note  here  that  Y-'  exists,  but  Q-1  does  not  exist  in  general.  It  is  at  this  stage  that  an 
important  point  should  be  noted;  in  general,  it  is  not  possible  to  reorganise  (29)  to 
achieve  T  and  W  on  one  side  of  the  equation  and  A  and  R  on  the  other  because  of  this 
inversion  roblem.  This  seems  to  be  assumed  possible  in  the  cascade  formulations  [4], [5] 
in  order  to  calculate  a  composite  transmission  matrix.  If  the  block  matrix  in  (29)  is 
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suitably  truncated  (ie  O  is  suitably  truncated)  it  may  be  possible  to  invert  Q  (and  thus 
reorganise  (29))  but  such  a  procedure  may  be  very  dangerous  because  the  original  O  is 
rank  deficient.  What  is  required  is  a  scheme  that  does  not  involve  the  inversion  of  Q. 
Cwik  and  Mittra  [3]  describe  one  such  method  while  the  authors  will  describe  another. 


BOUNDARY  CONDITIONS  AT  A  RESISTIVE  SHEET 

We  now  consider  the  second  canonical  interface  (2),  and  asume  a  resistive  sheet  of 
1/<7S  Ohms  per  square  between  two  electrical  media 
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In  this  case  there  is  no  cross  coupling  between  modal  fields  and  thus  a  surface  current 
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(31) 


is  impressed  on  the  surface  where 

E(t)  -  A  +  R  -  T  W  (32) 

pqr  pqr  pqr  pqr  pqr 

since  the  tangential  electric  field  is  still  continuous  across  the  interface. 

We  now  have  a  discontinuity  in  the  tangential  magnetic  field  given  by 
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thus,  assuming  the  previous  matrix  notation  where  {Ej('  .2)} 
and  (33)  may  be  combined  to  give  us 


{fpqr(1'2)}-  ‘he"  (31).  (32) 
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As  opposed  to  (29),  (34)  is  a  matrix  whose  elements  may  be  properly  inverted  if  required. 
However,  since  this  may  not  be  done  with  (29)  a  method  must  be  found  to  connect 
matrices  of  this  form  which  does  not  require  a  reorganisaiton.  Such  a  method  is 
illustrated  below. 

CONNECTION  OF  MATRICES  WITHIN  A  GENERAL  CASCADED  SYSTEM 

A  cascade  system  that  does  not  require  the  inversion  of  the  matrix  O  in  (30)  may 
be  constructed.  This  has  the  further  advantage  that  evanescent  modal  fields  and  fields 

within  a  lossy  dielectric  can  be  organised  in  such  a  way  that  they  always  appear  to  decay 
in  the  numerical  implementation  for  arbitrary  thicknesses  of  dielectric.  This  means  that 
the  number  of  evanescent  connections  between  interfaces  can  be  highly  overspecified 
without  causing  any  ill-conditioning  of  the  scheme.  The  method  proceeds  by  a  grouping 
of  vwo  contiguous  canonical  interfaces  separated  by  a  dielecuiclmagnetic  spacer  to  form  a 
single  composite  structure  whose  form  is  also  that  of  (29)  and  (34).  This  composite 

structure  is  then  taken  together  with  the  next  canonical  interface  and  the  method  repeated 

in  a  systematic  fashion  over  the  entire  multi-layer  structure.  The  grouping  of  two 

contiguous  structures  is  illustrated  below  in  Fig  3. 
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Figure  3 


11 


Suppose 


(35) 


(36) 


where  Xj j '  and  Xjj"  are  known  but  non-invertible  and  Wj  and  are  unknown 

quantities.  In  general,  if  the  material  between  interface  (1)  and  (2)  is  fr,/ir  separated  by 
a  distance  d,  then 


hi} 

and 

^pqC 


(37) 

(38) 

(39) 


(40) 


defined  so  that  TpqC  (independent  of  £)  is  defined  with  Im(7pqf)  ^  0.  Equations  (35), 
(36),  (37)  and  (38)  may  be  combined  as 


-  Xil*l  +  X12  H  X22  h]  +  X12  H  X21  H  h 


(41) 


We  may  thus  write 

I,  -  (I  -  L)'1  a 

where 

a  -  xjjA,  +  xj2  h  x'22  w2 

and 
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Thus,  even  though  L  *  may  not  exist,  (I  —  L)  ^  does  exist  and  is  well  conditoned. 
We  can  now  obtain  T_2  an^  Si  as 
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(45) 
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Z12  -  Xil 
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L)'1  X'n  H  X"2  +  Xj 2 

(47) 

Z21  “  X21 

+  X22 

H  X^j  H  (I  -  L)'1  X'u 

(48) 

Z22  “  (X22  H  X21 

H(  1  -  L)'1  Xj2  +  X^}h  X22 

(49) 

With  T2  and  Rj  obtained  from  the  first  two  interfaces,  the  next  interface  is  treated  using 
its  own  conneciton  matrix  and  that  just  generated  via  (45)  to  construct  a  new  conneciton 
matrix  uisng  the  same  method  as  outlined  above.  This  proceeds  until  there  are  no  more 
interfaces  in  which  case  we  have  a  final  incoming  W  =  0  and  T  and  R  are  the  final 
required  vectors  of  the  composite  structure.  The  method  differs  from  that  of  Cwik  and 
Mittra  [3]  in  two  important  respects.  Firstly,  the  matrix  H  never  appears  in  an  inverted 
form  so  exponential  terms  never  cause  numerical  problems  by  becoming  too  large  and 
secondly,  our  approach  uses  electric  fields  rather  than  power  waves. 

In  this  formulation  it  may  immediately  be  seen  that  interactions  between  any  pair  of 
interfaces  in  fig  3  are  coupled  through  the  propagation  equations  (32)  and  (38).  This 
means  that  the  sizes  of  the  Tj,  A2,  Wj  and  R2  may  be  limited  in  size  to  encompass 
only  those  evanescent  modes  whose  decay  (described  by  H)  is  non-negligible  over  the 
distance  d.  We  may  further  note  that  Aj ,  Rj  and  T2,  Wj,  may  both  be  of  different 
lengths  in  the  situation  where  distances  between  interfaces  vary  greatly  in  size  (so  that  the 
same  number  of  evanescent  modes  is  not  required  to  couple  each  interface).  In  fact,  this 
will  generally  be  required  since  usually  only  the  zero-order  modal  coefficients  will  be 
required  in  the  air  region  to  either  side  of  the  composite  FSS.  This  will  then  result  in 
the  matrices  X^',  Xjj"  and  Zjj  and  Z22  being  rectangular. 

One  possible  scheme  that  may  be  imposed  is  to  set  an  automatic  threshold  0  <  T  < 
1  and  to  limit  the  matrix  size  so  that  only  those  modes  p,q,r  are  taken  in  determining 
the  size  of  Xl.  &2'  3^1  and  R2  such  that  exp[+jypq,4]  <  T  from  (39), (40). 

TRANSMISSION  AND  REFLECTION  S-MATR1CES  AND  ENERGY  BALANCE 

Once  the  final  T  and  R  are  known  in  terms  of  the  initial  vector  A,  the  problem  is 
almost  solved.  Usually,  the  final  quantities  required  are  the  transmitted  and  reflected 
electric  fields  in  the  $  and  £  directions  given  an  incident  wave  of  unit  amplitude  and 
arbitrary  polarisation  incident  at  the  angle  (c>0,90).  With  such  a  full  polarimetric 
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description  it  is  convenient  to  describe  a  polarisation  coordinate  system  in  the  direction  of 
the  wave. 

Let  us  write  the  incident  electric  field  by 


£  sin  £  x  +  cos  f  v 

*0  r0  ■  v0 

6  -  cos  f  cos  8  x  +  sin  <p  cos  8  v  -  sin  8  z 

o  yo  o  o  o*  o 


(50) 


(51) 


the  reflected  wave  has  an  angle  of  reflection  Prefl  =  60,  fTtf\  =  T  +  v°0'  but  the  wave  is 
directed  away  rather  than  towards  the  surface.  Thus  define 


6f  ^refl  +  ^8  ^refl 


(52) 

(53) 


With  the  transmitted  field,  the  transmission  angles  #tran  =  *  -  80  v’tran  =  x  +  anci 
the  wave  is  directed  in  the  same  direction  as  the  incident  wave,  thus  define 
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(57) 


Erefl 


rll  r12 
r21  r22 


It  may  now  be  shown  that  if  the  final  matrix  equation  (45),  produced  after  all  the 
interfaces  have  been  included,  is  given  by 


2  2 

Til  712 

2  Z 

Z21  z22 


with  the  p  =  q  =  0  terms  of  Zjj  given  by 


;o) 

Id) 

zn 

Z12 

;<n 

Z21 

;o) 

22 

=  1  and 

2  inte 

1(2) 

1(2) 

211 

2  12 

;(2) 

Z21 

;(2) 

z22 

(for  r  =  1  and  2  interactions)  and  with  the  p  =  q  =  0  terms  of  Z21  given  by 


(for  r  =  1  and  2  interactions) 


C11 

*12 

;(i) 

11 

~(i)  0  l 

-Zj2  cos  8 

l21 

l22  . 

-z(l> 

21/cos  0 

;a> 

z22 

and 

rll 

r12 

;(2) 

Z11 

-Zj^cos  e 

,  r21 

r22  . 

;(2) 

21/cos  e 

,;(2) 

Z22 

furnishing  the 

desired  result.  If  the 

composite  FSS 

(58) 


(59) 


conservation  of  energy  implies 


2  2  2 
'll  +  l21  +  rll  +  r2 1 


-  1 


(60) 


and 


2  ,  2  2 
1 12  *22  +  r21  +  r22 


-  1 


(61) 
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for  arbitrary  angles  of  incidence  provided  that  no  higher  order  grating  lobes  are  excited. 
This  provides  a  useful  check  on  the  calculation.  Unfortunately,  although  a  necessary 
condition  (60)  and  (61)  are  not  a  sufficient  condition  and  it  is  not  appropriate  to  use  (60) 
and  (61)  to  attempt  to  provide  an  estimate  on,  for  example,  the  number  of  modes 
required  in  the  program.  When  losses  or  grating  lobes  are  present  we  have  the  strict 
inequalities. 


<  1 


<  1 


get  significant  energy  absorption. 

MEASUREMENTS  FACILITY 

In  order  to  confirm  the  method,  predictions  using  this  program  and  experimental 
measurements  of  the  transmission  coefficient  using  a  Vector  Network  Analyser  have  been 
made  for  a  selection  of  frequency  selective  surfaces  at  various  angles  of  incidence  and 
polarisation.  The  experimental  arrangement  is  illustrated  below  in  Fig  4. 


l‘nl 
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+ 
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1 2  | 
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1  r  2 1  I 

2  | 
+ 

r22  1 

and  with  lossy  dielectrics  it  is  often  possible  to 


16 


Measurements  have  been  successfully  made  from  4.5  GHz  to  20  GHz  in  ranges 
determined  by  the  horns  used  to  transmit  and  receive.  The  FSS  is  about  4x4  feet  in 
size,  rotatable  about  its  vertical  axis  and  the  horns  may  be  rotated  so  that  the  E -field  is 
either  vertical  or  horizontal.  The  separation  distances  Lj  and  L2  are  of  the  order  of  8 
feet  providing  an  approximate  far-field  plane  wave  illumination  of  the  screen. 
Measurements  are  made  with  and  without  the  FSS  panel  and  the  appropriate  vector 
division  is  carried  out  automatically  by  the  VNA  providing  both  power  transmission  and 
phase  (if  required)  coefficients  of  the  FSS. 

PREDICTIONS  AND  EXPERIMENTAL  COMPARISONS 

In  the  following  examples  rectangular  slotted  structures  are  considered.  The  reasons 
for  this  are  that  (1)  their  modes  and  thus  the  best  choice  of  field  expansion  modes 
jipqr(x)  are  well  known  (see,  for  example  Chen  [6])  and  (2)  they  are  of  considerable 
practical  interest  for  the  design  of  frequency  filters  that  are  also  required  to  be 
polarisation  isolating. 

Example  1 

Consider  a  single  array  of  slots  backed  onto  a  dielectric  substrate.  The  array  of 

slots  has  a  unit  cell  size  (see  fig  1.1)  with  dx  =  17.5mm,  dy  =  7.5mm,  o  =  41 '  with  a 
single  slot  per  unit  cell  of  length  15.0mm  by  1.0mm  aligned  with  respect  to  the  unit  cell 
so  that  its  long  dimension  is  pointing  in  the  x-dirn.  The  dielectric  backing  is  standard 
fibre-glass  FR4  material  of  thickness  1.57mm  with  an  estimated  Re(er)  =  4.0  and  loss 

tangent  tan  &  =  0.027  (These  values  are  assumed  in  all  the  numerical  predictions  by  the 
authors  using  FR4,  and  seem  the  best  values  to  take  over  the  frequency  range  7.0  GHz  -> 
17.0  GHz). 

In  the  following  predictions,  the  number  of  Floquet  modes  used  are  given  by  Pmax. 

=  Omax  =  13,  and  the  expansion  modes  ^pqr  are  chosen  as  ^pqr  =  gp^  for  q  =  0  and 

-10  *  p  «  10.  (This  is  not  a  good  choice  for  j^pqr  with  slots,  but  is  sufficient  to 
illustrate  the  method).  Finally,  the  number  of  evanescent  modes  required  to  couple  the 
air-dielectric  and  slotted  array  interface  is  taken  with  an  automatic  threshold  of  T  =  0.05 
at  all  frequencies  and  angles  of  incidence.  Figures  5.1  a-c  illustrate  the  predictions  and 
measurements  of  the  power  transmission  coefficient  in  db,  given  by  201og]g|tjj  1,  over  the 
range  7  to  17  GHz  for  transverse  electric  polarisation  (E  directed  in  the  y-direction)  for 
6  =0,  10  and  45  ,  (y>  =  0).  The  continuous  curves  illustrate  the  measurements  and  the 
dots  the  predictions.  Figures  5.1  d-e  illustrate  the  predictions  and  measurements  for  (*  = 
10’,  45’  (y>  =  90')  for  the  transverse  magnetic  polarisation  (H  directed  in  the  x-direction) 
given  by  20  •og1git22l. 

Example  2 

Consider  the  same  array  of  slots  as  used  above,  but  with  one  on  either  side  of  the 
dielectric  substrate.  The  two  arrays  are  arranged  so  that  there  is  a  lateral  offset  between 
them  in  the  ^-direction,  but  no  offset  in  the  y-direction.  If  100%  offset  indicates  a 
lateral  shift  in  a  given  direction  by  a  complete  unit  cell,  then  in  this  example  there  is  a 
50%  lateral  offset  of  the  slots  in  the  x-direction  (along  the  axis  of  the  slot). 

Figures  5.2  a-c  illustrate  the  predictions  and  measurements  for  transverse  electric 
polarisation  at  6  =  0,  10  and  45’,  and  figures  5.2  d-e  illustrate  the  predictions  and 
measurements  for  transverse  magnetic  polarisation  at  f>  =10  and  45'  as  in  example  1. 
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Example  3 


Consider  3  arrays  of  slots,  of  the  kind  described  above,  separated  by  two  thicknesses 
of  FR4  material  both  of  width  1.57mm.  In  this  example,  all  3  slot  arrays  are  aligned 
with  each  other  with  no  lateral  offset.  The  frequency  range  examined  is  from  4.5  GHz 
to  12  GHz  and  a  figure  for  the  threshold  T  is  chosen  as  before  with  T  =  0.05  (where  it 
may  be  noted  that  about  41  transverse  electric  and  41  transverse  magnetic  connecting 
modes  are  necessary  at  12  GHz).  Figures  5.3  a,b,c  show  predictions  and  measurements  at 
6=0,8  =  10"  TE  and  6  =  10'  TM  in  the  sense  defined  in  example  1. 

CONCLUSIONS 

The  method  described  is  the  most  robust  multi-layer  FSS  method  that  the  authors 
are  aware  of  with  significant  accuracy  advantages  over  some  other  cascade  methods. 
Provided  a  sufficient  number  of  Floquet  modes  is  chosen  together  with  a  sufficiently  good 
set  of  field  expansion  functions  (necessary  criteria  for  any  single  layer  FSS  program)  then 
the  method  described  is  well  conditioned  and  capable  of  high  accuracy  when  (1)  there  are 
significantly  fewer  field  expansion  functions  than  evanescent  modes  and  (2)  the  coupling 
threshold  T  is  chosen  arbitrarily  small.  The  examples  presented  are  by  no  means  special 
and  probably  do  not  have  any  special  practical  application,  but  a  more  detailed  parameter 
study  is  currently  underway  to  examine  slotted  arrays  of  more  interest.  This  will  be 
reported  on  in  a  later  article. 
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Figure  5.1a  FSS  Examp!e_1  theta=0  degrees  .CTEMD 
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Figure  5.1b  FSS  Example_1  theta=10  degrees  ,CTE) 
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Figure  5.1c  F5S  Exampie_1  theta=45  degrees  ,CTE) 


Figure  5. Id  FSS  Examp!e_1  theta=10  degrees  ,(TM) 
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Figure  5.2b  FSS  Example_2  theta=1D  degrees  ,[TE3 
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Figure  5.2c  FSS  Example_2  theta=45  degrees  ,(TE) 
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Figure  5.3a  FSS  Example_3  theta=0  degrees  ,(TEM) 
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Figure  5.3c  FSS  Exampie_3  theta=10  degrees  ,(TM] 
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